global par_;

switch switch_distribution
    case 'gaussian'
        [T2, Pgrid] = TRP_AR(rho, 0, par_.std_eA^2, nump, -0.5, 0.5);
    case 'laplace'
        [T2, Pgrid] = TRP_AR_laplace(rho, 0, par_.std_eA^2, nump, -1, 1);
    case 'mixed'
        [T2, Pgrid] = TRP_AR_mixed(rho, 0, par_.std_eA^2, 0.9, lam, nump, -1, 1);
end

PMAT = exp(Pgrid);  % price grid in levels

pstep = Pgrid(2)-Pgrid(1);  % distance between price grid points

R = Rmatrix(nump, mu, pstep);
